Exactly solvable case of a one-dimensional Bose-Fermi mixture 



in 
o 
o 

> 
O 



i 

CO 

c3 



I 

C 

o 
o 



(N 
> 

m 

in 
o 
in 
o 

03 



S3 
O 

o 



X 
S3 



Adilet Imambekov and Eugene Demler 
Department of Physics, Harvard University, Cambridge MA 02138 
(Dated: February 2, 2008) 

We consider a one dimensional interacting bose-fermi mixture with equal masses of bosons and 
fermions, and with equal and repulsive interactions between bose-fermi and bose-bose particles. 
Such a system can be realized in experiments with ultracold boson and fermion isotopes in optical 
lattices. We use the Bethe-ansatz technique to find the ground state energy at zero temperature 
for any value of interaction strength and density ratio between bosons and fermions. We prove that 
the mixture is always stable against demixing. Combining exact solution with the local density 
approximation we calculate density profiles and collective oscillation modes in a harmonic trap. 
In the strongly interating regime we use exact wavefunctions to calculate correlation functions for 
bosons and fermions under periodic boundary conditions. 



Recent developments in the cooling and trapping of 
atomic gases open exciting opportunities for experimen- 
tal studies of interacting systems under well-controlled 
conditions. Using Feshbach resonances and/or opti- 
cal lattices 0, El it is possible to reach strongly inter- 
acting regimes, where correlations between atoms play 
a crucial role. The effect of interactions is most promi- 
nent for low dimensional systems, and recent experi- 
mental realization 0, Q of a strongly interacting Tonks- 
Girardeau (TG) gas of bosons opens new perspectives 
in experimental studies of strongly interacting systems 
in ID 6]. A different line of research, which has at- 
tracted considerable attention both theoretically^ and 
experimentally 8, 9], is the study of atomic mixtures. ID 
Fermi gases have been theoretically well investigated [l"oj. 
since they correspond to interacting electrons in solid 
state quasi-lD systems. On the other hand, ID bose- 
fermi(BF) mixtures did not attract equal attention until 
recently Il2t Il3t Il4j . when it became possible to real- 
ize such systems in experiments with cold atoms. Several 
properties of such systems have been investigated so far, 
including phase separation, fermion pairing and charge 
density wave (CDW) formation. However, most of these 
investigations relied on the mean-field approximations or 
the Luttinger liquid (LL) formalism. The mean-field ap- 
proximation is known to be unreliable in ID, and the LL 
approach does not allow the calculation of LL parame- 
ters in the strongly interacting regime. In this paper we 
study BF mixtures in the regime where an exact Bethe- 
ansatz solution is available. We use the exact solution to 
calculate the ground state energy and investigate phase 
separation and collective modes. In the strongly interact- 
ing regime we derive an analytical determinant formula, 
which allows to calculate correlation functions at all dis- 
tances numerically for a polynomial time in the system 
size. The exactly solvable case considered in this paper 
is relevant to current experiments, and can be used as a 
benchmark to check the validity of different approaches. 



A ID interacting BF mixture is described by the 



Hamiltonian 



H= f dx{^—d x ^>\d x ^ b + ^—dMd x ^ s ) + 
Jo 2m b 2m/ s 

J ds(-flM¥j¥j*6** + S&/*Jtf >*/*&)• (!) 

Here, ^b,^f are boson and fermion operators, mi,, rut 
are the masses, and gbb-, gbf are bose-bose and bose-fermi 
interaction strengths. An array of such one dimensional 
tubes can be realized using a strong optical lattice in y, z 

directions!! HE!- 

Alternatively, one can create a single 
BEC in a ID box potential using the techniques of Ref. 
(lif. The model (JIJ is exactly solvable, when 



m f = m b = m, g bb = gbf = g. 



(2) 



The first condition is approximately satisfied for iso- 
topes of atoms. Some of the promising candidates are 
m Y b+ 172 Yb^ z^K-^K,!?] and 86(84)^ _87(M) 
Kb 19] . Alternatively, if one uses an additional op- 
tical lattice along the x direction with filling factors 
smaller than one, then (JIJ is an effective Hamiltonian de- 
scribing this system with the masses determined by the 
tunneling jfj|. Thus one can simulate Q using already 
available degenerate mixtures 8]. We point out that to 
satisfy the second condition in (J2J, it is sufficient to 
have equal (positive) signs for the two scattering lengths, 
but not necessarily their magnitudes. Well away from 
confinement induced resonances (2^, ID interactions are 
given by g b b = 2haj b ±a b b, gbf = 2h^/Lu b ±ujfj_abf, where 
LUb±, u>f± are radial confinement frequencies, and abb, a bf 
are 3D scattering lengths. For a fixed value of a b b/ob b f, 
one can always choose the detuning of the optical lattice 
laser frequencies in such a way that gbb — gbf- Addi- 
tionally, all interactions can be tuned using available BF 
Feshbach resonances 

Model under conditions (J2J has been considered in 
the literature before [2l|, but its properties have not been 
investigated in detail. Similar to the case of bosons , 
a dimcnsionless parameter which controls the strength of 
interactions is 7 = mg/(h 2 n), with n being the total den- 



2 



sity of particles. 7 > 1 corresponds to the strongly in- 
teracting mixture. The ground state energy of the model 
||TJ, J2J with periodic boundary conditions is related to 
the solution of a system of integral equations 



27T /9 (fc) = 1 + / 

2na{A) -- 



B 4ca(A)dA 
B c 2 + 4(A - fc) 2 ' 
Q Acp(uj)duj 
Q c 2 +4(A- w) 2 ' 



(3) 
(4) 



Here c = mg/H 2 , and Q, B are related to densities via 

rQ i-b 
N/L = / p(k)dk, N b /L = / cr(A)dA, (5) 

J-Q J-B 

where Nb is the number of bosons, N is the total num- 
ber of particles, and L is the length of the system. The 
ground state energy can be expressed as 



h 2 r Q 

E/L = — / k 2 p(k)dk. 
2m J-q 



(0) 



One can solve eqs. ll.')l- r >l) numerically and calculate the 
energy as a function of 7 and relative boson density 
a = N b /N. After rescaling E = e(7, a)h 2 N 3 / (2mL 2 ) , 
where e(7,a) is shown in fig^ I n the limit 7 > 1, one 
can derive the large 7 expansion of the energy e(7, a) = 
(7r7) 2 /3(7 + 2(a-|-sin(7ra)/7r)) 2 . Chemical potentials for 
bosons and fermions are given by pb = 9E / dNb, fif — 
8E/dNf, where Nf — N — Nb- Using the exact solu- 
tion, one can analyze demixing instabilities [ill ll2T | for 
repulsive BF mixtures. In the absence of an external 
potential the BF mixture is stable if the compressibility 
matrix dpb(f)/d n b(f) is positively defined(here, n.&(/) is a 
boson(fcrmion) density). In our system we verified nu- 
merically that this condition is satisfied for all a and 7. 
This proves that the BF mixture under conditions J5J is 
stable against demixing for any strength of interaction. 
This is in contrast to mean field approximation which 
predicts demixing for strong enough interactions |l l| . Al- 
though an exact solution is available only under condi- 
tions (J2J) , small deviations from these should not dramat- 
ically change the energy e(7, a). Therefore, we expect the 
ID mixtures to remain stable to demixing in the vicinity 
of the integrable line (J2J for any interaction strength. 



Using numerically obtained e(j, a) and local den- 
sity approximation(LDA), one can investigate the be- 
havior of the system in an external field. Local den- 
sity approximation is expected to be valid for slowly 
varying external potential. We will assume that the 
confinement frequency for bosons and fermions is the 
same. Within LDA, density distributions in the re- 
gion of coexistence are governed by a set of equa- 
tions fi f {n° b (x),n° f (x)) + = /^(n°(0),n°(0)) and 



p b (n b (x),n° f (x)) + ': 



fif (ng (0) , n° f (0)) . It turns out 



that this set of equations cannot be satisfied for the whole 
cloud, and consequently a phase separation occurs in the 
presence of such external potential. The BF mixture is 
present in the central part, but the outer sections consist 
of Fermi gas only. In the weakly interacting limit, this 
can be interpreted as an effect of the Fermi pressure 0. 
While bosons can condense to the center of the trap, 
Pauli principle pushes fermions apart. As interactions get 
stronger, the relative distribution of bosons and fermions 
changes, and fig|21contrasts the limits of strong and weak 
interactions. For strong interactions, the fermi density 
shows strong non-monotonous behavior. 

Recent experiments^^ demonstrated that collective 
oscillations of ID gases provide useful information about 
interactions in the system. Within the region of coex- 
istence of bosons and fermions, such oscillations can be 
described by four hydrodynamic equations |23J 



m 



d d 1 

'dt Svb{f} + dx~^ Hf} + Vext ' b ^ + 2 TOUfc2 (^ = 



d 



0, (7) 
0. (8) 



In certain cases, analytical solutions of hydrodynamic 
eqs. are available |23l |24| and provide the frequen- 
cies of collective modes. When an analytic solution 

is not available^ the "sum rule" approach has been 

used[H |H |H |2f|. The disadvantage of the latter ap- 
proach is an ambiguity in the choice of the multipole 
operator which excites a particular mode, especially for 
multicomponent systems 26J. Here we develop an effi- 
cient numerical procedure for solving the hydrodynam- 
ical equations in ID, which doesn't involve additional 
"sum rule" approximation. After linearization, a system 
of hydrodynamic equations can be written as 



Sfj.b(x) 
Sfj, f (x) 



dn b 
dn b 



dp.b 
dnj 

dn f 



n° b (x) 
n°Jx) 



r 



5pb{x) 
Sfj, f (x) 



(9) 



In the outer fermionic shell, <5/ij"* satisfies an equation: 



mu^iif* = ^-V \n° f ut (x)VSp° f ut } 



(10) 



All modes can be classified by their parity with respect 
to the space inversion, x — > — x. Thus, we need to solve 
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equations (|9I10() only in the region x > 0. For even 
modes, one has the following conditions at x = : 

VS/jif(x = 0) = 0, V5m&0 = 0) = 0. (11) 

For odd modes, analogous conditions are 

5fj, f (x = 0) = 0, 6fi b (x = 0) = 0. (12) 

Boundary conditions for fermions at the edge of the 
bosonic cloud, x b , correspond to the continuity of Sfif 
and Vf : 

VSfif^x = Xb + 0)= V8nf(x = x b -Q), (13) 
6(i° f ut (x = x b + 0) =Sfj,f(x = x b -0). (14) 

Two additional conditions come from the absence of 
the bosonic(fermionic) flow at x b (xf) : 

nl{x)v b {x)\ x ^ Xb - = 0, (15) 
rif(x)vf(x)\ x -> Xf -o = 0. (16) 

Outside of the region of coexistence, the chemical po- 
tential and density of fermions are given by /iy nt ~ 

(n™*) 2 ,riy nt ~ y/1 — (x/xf) 2 , where Xf is the fermionic 
cloud size. For these functions, there exists a gen- 
eral nonzero solution of eq. I|10[l which satisfies (|16|l : 
Sfi°j ut = cos((o;/wo) a,rccos x / x f) . 

Substituting this into l)13 [l -l)14 [l . one has to solve eigen- 
mode equations numerically for x < x b , with five bound- 
ary conditions |JT3J| , {T1J , and {TTJ or {115} depending 
on the parity. These boundary conditions are compati- 
ble, only if u> is an eigenfrequency. We choose to leave out 
condition (|15|) . find a numerical solution to a system of 
differential equations, and check later if (|15fl is satisfied 
to identify the eigenfrequencies. 

First we apply this numerical procedure for weakly- 
interacting regime, see the inset of fig. |3 When 70 — > 0, 
bose and fermi clouds do not interact, and collective 
modes coincide with purely bosonic or fermionic modes, 
with frequencies [23j = nujQ and tu b = ujQ^/n(n + l)/2. 
As interactions get stronger, bose and fermi clouds get 
coupled, and all the modes except for dipole Kohn mode 
change their frequency. For the case of strongly interact- 
ing BF mixture, we find low-energy modes which corre- 
spond to the "out of phase" oscillations of bose and fermi 
atoms, and that keep the total density approximately 
constant. These modes can be understood as follows: for 
7 ^> 1 the energetic penalty for changing the relative 
density of bosons and fermions is small. The frequency 
of such relative modes scales as ~ wo/^/to- In fig. 01 we 
show the numerically obtained dependence of the energy 
of these low energy modes as a function of the ratio of 
the total number of bosons and fermions. 

Finally, we report the results of the calculations of 
the bose-bose and fermi-fermi correlation functions in 
systems with periodic boundary conditions. Such cal- 
culations can be performed in the regime 7 » 1, due 




FIG. 1: Energy of the ground state is given by E = 
e(y,a)h 2 N 3 /(2mL 2 ), where 7 = mg/(h 2 n), and a = N b /N is 
the boson fraction. When a — 0, system is purely fermionic, 
and the energy doesn't depend on interactions. When a = 1, 
the system is purely bosonic, and numerically obtained en- 
ergy coincides with the result of [2^| . If 7 = 0, bosons and 
fermions don't interact, and e(j,a) — (7r 2 /3)(l — a) 3 . 

to "factorization" of the Bethe wavefunction, similar to 
[2?|. "Spin" part of the wavefunction is related to the 
ground state of spin XY model, and significant analytical 
progress can be made compared to the case of spin— 1/2 
fermion model |2jj . Details of these calculations will be 
presented elsewhere |28j. The Fourier transform of cor- 
relation function is an occupation number n(k), which 
can be measured directly in time-of-flight experiments 
or using Bragg spectroscopy 29|. For fermions, n(k) is 
shown in fig 01 The discontinuity at kf — irnj gets 
smeared out by interactions, but the overall change in 
n(k) after crossing kf remains quite large. For bosons, 
n(k) is a monotonous decreasing function, with a singu- 
larity at k = 0. Singularity at k = is governed by the 
long distance behavior of correlation function, character- 
ized by the bosonic LL [33,1101 parameter K b : 

p{0,O-l/\e^ Kb) Mk)r.\k\-'+y^\ (17) 

where £ is a distance between two points. For periodic 
boundary conditions, one can extract K b using [3 lj a fit 
p(0,0 - l/|sin(7rC/i)| 1/(2iiri>) for £ > L/N b . Results 
of the numerically extracted value of K b for strong in- 
teractions are shown in fig ^ depending on the rela- 
tive number of bosons and fermions, K b changes from 
0.5 to 1. We note that in the absence of fermions K b 
changes from 00 to 1 for interacting bosons as the inter- 
actions get stronger. In the presence of fermions, bosonic 
correlations decay faster than in the purely bosonic sys- 
tem. In principle, the low energy physics of a BF mix- 
ture is characterized not only by bose-bose or fermi-fermi 
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FIG. 2: Densities of Bose and Fermi gases in weakly interact- 
ing(left , 70 = 0.18) and strongly interacting^ S> 1) regimes. 
The overall number of bosons equals the number of fermions. 
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FIG. 3: Dependence of the energy of several lowest lying 
"out of phase" modes for 70 2> 1 on overall boson fraction a, 
where 70 is the TG parameter in the center of a trap. Note, 
that we don't show "in phase" TG modes uj = nuo, since 
their frequency is much higher. The inset shows frequencies 
of low-lying modes in mean-field regime for an equal overall 
number of bosons and fermions. 



correlations, but also by the long-range behavior of BF 
correlations [l^ LL-I- 0™ method can be applied to the 
evaluation of BF correlations and can be generalized to 
nonzero temperatures. 
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